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ABSTRACT 

Recent data on Galactic cosmic-ray (CR) leptons and hadrons gave rise to two exciting problems: on the 
lepton side, the origin of the rise of the CR positron fraction e+/(e“ -i- e"*") at ^ 10 - 300 GeV of energy; on the 
hadron side, the nature of the spectral hardening observed in CR protons and nuclei at ^ TeV energies. The 
lepton anomaly indicates the existence of a nearby source. It has been proposed that high-energy positrons 
can be produced inside nearby supernova remnants (SNRs) via interactions of CR hadrons with the ambient 
medium. A distinctive prediction of this mechanism is a high-energy rise of the boron-to-carbon ratio, which 
has not been observed. It also requires old SNRs at work (with ineffective magnetic field amplification and 
slow shock speed), that cannot account for the CR hadronic spectra observed up to the knee energies 5 PeV). 
We propose a new picture where, in addition to such a nearby CR accelerator, the high-energy spectrum of CR 
hadrons is provided by the large-scale population of SNRs, on average younger, that can efficiently accelerate 
CRs up to the knee. Under this scenario, the spectral hardening of CR hadrons can be naturally interpreted as 
the transition between the two components. As we will show, our two-component model breaks the connection 
between the positron fraction and the boron-to-carbon ratio, which is now predicted to decrease with energy in 
accordance with the data. Eorthcoming data from AMS will be crucial for testing this model. 

Subject headings: cosmic rays — acceleration of particles — ISM: supernova remnants — nuclear reactions, 
nucleosynthesis, abundances 

1. INTRODUCTION 

The AMS collaboration has recently confirmed with high 
precision the CR positron fraction anomaly previously ob¬ 
served by PAMELA (Adriani et al. 2009; Aguilar et al. 2013). 

The data show a rise of the fraction up at energies ~ 10 - 
300 GeV, followed by a possible plateau at higher energies 
(Accardo et al. 2014), which cannot be described by conven¬ 
tional models of e+ production by collisions of CR hadrons 
with the interstellar matter (ISM). In these models, the pri¬ 
mary CRs {e.g., e~ , H, He, C, or O) are accelerated by 
supernova remnants (SNRs) via diffusive shock acceleration 
(DSA) mechanisms up to ^ PeV energies to power-law spec¬ 
tra S cx E~'^. Their subsequent propagation is described by 
means of an energy-dependent confinement time t^sc oc E~^ 

(or diffusion coefficient K oc E^), and their collisions with 
the ISM give rise to secondary CRs such as Li, Be, B, e+, 
or p (Strong et al. 2007). At E 2^ GeV, this picture pre¬ 
dicts power-law spectra for primary nuclei, Np ^ E~'^~^, 
and decreasing secondary-to-primary ratios, e.g., the boron- 
to-carbon: B/C ^ E~^, where <5 ^ 0.3-0.7, u ~ 2.0-2.4, and 
1 / + S 2.7. The high-energy spectrum of CR leptons is fur¬ 
ther steepened by radiative losses, with characteristic time- 
scale t’^‘^‘^{E) oc E~^. Given all the positrons are from sec¬ 
ondary origin, the positron fraction is expected to decrease 
similarly to other secondary-to-primary ratios, in remarkable 
contrast with the observations. Explanations of the rise may 
include either dark-matter particles annihilation or decay, or 
nearby astrophysical sources (Serpico 2011). Among the sec¬ 
ond class, it has been proposed that high energy positrons 
may be produced through hadronic interactions of CR protons 
undergoing DSA inside old SNRs (Blasi 2009). Yet, if sec¬ 
ondary positrons are produced and accelerated by this mech¬ 
anism, other secondary species (such as B nuclei or p) will 


also be produced from CR nuclei interacting with the gas. As 
shown by Mertsch & Sarkar (2009), this mechanism leads to 
an increase of the B/C ratio at > 100 GeV per nucleon. How¬ 
ever the current B/C ratio data decrease with energy, indicat¬ 
ing that the old SNR scenario (hereafter OSNR) should be 
ruled out (Cholis & Hooper 2014). 

On the other hand, the spectra of primary CR elements have 
been measured up to ~PeV energies and beyond. Recent 
data on CR protons and nuclei revealed a remarkable spec¬ 
tral hardening at ^ TeV energies which stimulated great in¬ 
terest (Blasi 2013). According to the PAMELA data on H 
and He (Adriani et al. 2011), the change of slope is located 
at R Ri230GV of rigidity {i.e., momentum-to-charge ratio) 
with a very sharp transition, which is not seen by other exper¬ 
iments. While the sharpness of such breaks is under debate, 
the CR spectral hardening at > 1 TeV per nucleon is es¬ 
tablished by several experiments such as CREAM and ATIC 
(Panov etal. 2009; Yoon et al. 2010). The proposed expla¬ 
nations invoke acceleration mechanisms (Ptuskin et al. 2013), 
diffusion effects (Tomassetti 2012; Blasi, et al. 2012), or 
the superposition of local and distant sources (Bernard et al. 
2013; Thoudam & Horandel 2013; Vladimirov et al. 2012). 

In this Letter, we argue that the OSNR scenario is incom¬ 
plete in order to account for the observations of CR hadron 
spectra. Whether the OSNR represents a single source or 
a population of sources with the characteristics required for 
producing secondary {i.e., low shock speed and damped 
magnetic fields), it cannot provide the flux of CR hadrons ob¬ 
served in the ^ TeV - PeV energy region, so that an addi¬ 
tional high-energy component of CR accelerators is needed. 
We propose a two-component SNR scenario where the high- 
energy part of the CR flux is provided by a Galactic ensem¬ 
ble of SNRs, hereafter GSNR, that are on average younger 
and more efficient when accelerating primary hadrons at high 
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energies (but are unable to accelerate secondaries). A key 
consideration is that the local flux of light CR nuclei de¬ 
pends on the large-scale structure of the Galaxy, reflecting 
the contribution of a large population of SNRs and their his¬ 
tories (Taillet & Maurin 2003). For E >10GeV, their es¬ 
cape rate, Tg”^, is generally larger than their spallation rate 
in the ISM so that their propagation is limited only by 
the size of the diffusion region. In contrast, light leptons are 
subjected to radiative cooling, due to synchrotron radiation 
and inverse Compton scattering, which limits significantly the 
range that they can travel before reaching the solar system 
(Delahaye et al. 2010). Their characteristic distance is ap¬ 
proximately A ~ Vr’'®-diC oc confined to ~ 1 kpc at 

energy above ^ lOOGeV (depending on the propagation pa¬ 
rameters), which would legitimate the OSNR approach. Thus, 
while the observed can be largely produced by only one 
or few nearby sources, the total spectrum of CR protons and 
nuclei may well result as the sum two SNR components: the 
nearby OSNR component, which would dominate the flux be¬ 
low ~ lOOGeV, and the GSNR ensemble, which would pro¬ 
vide the high-energy flux up to the knee. As we will show, 
a two-component model gives a good description of the pri¬ 
mary CR spectral hardening and it has an impact on the spec¬ 
tral shape of the B/C ratio, which is now determined as the 
superposition of several contributions. 


2. CALCULATIONS 

The spectrum of CRs accelerated in SNRs is computed 
analytically using the linear DSA theory and including the 
secondary production terms due to hadronic interactions. 
We follow our calculation scheme in Tomassetti & Donato 
(2012), which has been proven to be equivalent to that of 
Mertsch & Sarkar (2009, 2014). In the shock rest-frame 
(x = 0), the upstream plasma flows in from a: < 0 with 
speed Ml (density ni) and the downstream plasma flows out 
to x > 0 with speed U 2 (density 712 ). The compression ratio is 
r = uilu 2 = 712 / 711 . For a nucleus with charge Z and mass 
number A, the DSA equation reads: 


9f Idu df ^ 

""dx ~^dx^ ^ 3 dx^ dp 


inel 


/ + Q j 


( 1 ) 


where / is the phase space density, D{p) is the diffusion 
coefficient near the SNR shock, u is the fluid velocity and 
rmei _ destruction rate for fragmenta¬ 

tion with the interaction cross-section The ambient den¬ 
sity n is assumed to be composed by 90% H and 10% He, 
similarly to the average ISM, in both sides of the shock. The 
source term is represented by Q. For primary nuclei Q'’*' = 
Y5{x)5{p — i.e., the ambient particle injection occurs 

immediately upstream the shock at momentum = ZR™\ 
where i?™-’ = 1GV for all species. The F-constants are abun¬ 
dance factors, determined from the data. The source term for 
secondary fragments produced by spallation of heavier (fc- 
labeled) nuclei is of the type where the 

partial rates of fragmentation are = cjdna^^^^. Using 
the subscript 7 = 1 (7 = 2) to indicate the quantities in the up¬ 
stream (downstream) region, the downstream solution of each 
nucleus can be expressed as: 


/ 2 (x,p) = /o + -(Qo-rr'/o) , (2) 

U2 


where the subscript 7 = 0 indicates the quantities evaluated at 
the shock (x = 0), and fo(p) is given by: 


Mp) = a r (P^y QP^^(p')e-xipY)^. 
Jo \PJ P 


+a 




(3) 


P' 


with a = 3r /(r — 1) and x ~ a(r™®yr’^®®)[i9(p) — D{p')], 
where F®^®® is the acceleration rate. The first term of Eq. 3 
gives the acceleration spectrum of primary particles at the 
shock, of the form ~ p~°‘e~^. The second term de¬ 
scribes the production and acceleration of CR fragments and 
it is coupled with the equations of heavier nuclei. The amount 
of secondary nuclei production depends on the SNR prop¬ 
erties via niDu^^. Their spectrum is of the form /q®® ~ 
/g”Z)(p)e“^ oc i.e., D times harder than the 

primary source spectrum. For having an efficient acceler¬ 
ation of all particles, the condition F'"®* <C F®^®® must be 
fulfilled in the whole momentum (or rigidity) range con¬ 
sidered (Mertsch & Sarkar 2009). For Bohm-type diffusion 
(D = ^), this condition can be also expressed as i? ^ i?® ^ 

20 ■ At i? > i?® one has F™®^ > F®®®, i.e., destructive 
interactions dominate over acceleration. In case of no inter¬ 
actions (<5 = 0 and r‘®®'Arag _ Q) ]-}jg usual DSA solution 
/ oc p~°‘ is recovered for primary nuclei, while the secondary 
CR production vanishes. The total CR flux produced by the 
SNR is obtained by integration in its volume: 




47rp^e 


max 



Attx^ f 2 j{x,p)dx 


(4) 


The exponential cut-off is used to explicitly account 

for the maximum momentum attained by the SNR, due to the 
finite time of the DSA process and it is assumed to oc¬ 
cur at the same rigidity for all CRs in the accelerator: = 

pinax^^ j^max |^g j-Qughly estimated from equating F®®® 

with 1/r™'', which gives R’^^^ 0.2i3Mfc“^r™'', though 

it is usually left as free parameter determined from the data 
(Kachelriefl et al. 2011; Serpico 2011). We stress that the 
steady-state description given here is an effective simplifica¬ 
tion of a more complex physical problem where the shock 
properties evolve with time. 

For the Galactic propagation, we use an analytical approach 
of CR diffusion and nuclear interactions, where the effects 
of energy changes are disregarded (Maurin et al. 2001). The 
Galaxy is modeled as a disc of half-thickness h containing the 
ISM gas (number density riism) and the CR sources. The disc 
is surrounded by a cylindric diffusive halo of half-thickness L 
and zero matter density. For each CR nucleus, the transport 
equation reads: 


dN _ d^N 
dt dz^ 


2hS{z)r’^^^N + 2h6{z)S, 


(5) 


where N{z) is its number density as function of the z- 
coordinate, K is the Galactic diffusion coefficient and r“®' = 
/ScTiismO"'®®* is the destruction rate in the ISM at velocity 
/3c and cross section cr™®/ The source term S includes the 
SNR contributions, and term for secondary production 
in the ISM from spallation of heavier (k) nuclei: 5'®™ = 
The diffusion coefficient is taken as K{R) = 






Connection between e“''/(e“'' + e ) and CR Hadron Spectra 


3 



Fig. 1.— Energy spectrum of CR positrons multiplied by E ^. The three 
models of Mertsch & Sarkar (2014) (lines) are compared with the new data 
from AMS (Aguilar et al. 2014). 

(3Ko{R/RoY, spatially homogeneous, where Kq expresses 
its normalization at i?o =4 GV. We solve Eq. 5 for all nu¬ 
clei (from Fe to FI) after assuming stationarity {dN/dt = 0), 
boundary conditions (N{zLL) = 0), and continuity condi¬ 
tion across the disc. The differential fluxes at Earth are given 
by (j){E) = ^No{E), where Nq, evaluated at z = 0, is of 

the type Nq ^ S . The quantities N, K, S 

and r depend on energy or rigidity too. To account for the 
solar modulation, we employ the force-field approximation 
(Gleeson & Axford 1968) using the parameter $ = 500 MV 
for a medium-level modulation strength. 

3. RESULTS AND DISCUSSION 

There are many parameters that determine the OSNR 
source spectra. We follow the benchmark model of 
Mertsch & Sarkar (2014), that provides good fits to the AMS 
leptonic data, assuming that the bulk of the flux is pro¬ 
duced by this type of OSNRs. All relevant parame- 

TABLE1 

Source and transport parameter sets . 


OSNR parameters Propagation parameters 


UI 

SxlO'^cms ^ 

Ko 

0.1 kpc^ Myr“^ 

B t Kb 

l/iG/16 

<5 

0.50 


4.65/4.55 

L 

5kpc 

n\ 

2cm-3 

h 

0.1 kpc 

^max 

ITV 

'^ism 

1 cm“® 

^snr 

50kyr 

<!> 

0.5 GV 


ters are listed in Tab. 1 . In particular we adopt B = l fj.G, 
i?max = 1 XV, Kb =16, and ui =5x 10^cms“^, where 
Kb parametrizes the deviation of D{p) from the Bohm 
value due to magnetic damping. These values are typi¬ 
cal for SNRs at their late evolutionary stages. The au¬ 
thors in Mertsch & Sarkar (2014) considered also scenar¬ 
ios with higher values of fixed at 3 TV and 10 TV, 

which can in principle discriminated with e’*' data at higher 
energy. In Eig. 1 we compare these predictions with the 
new high-energy data released by AMS (Accardo et al. 2014; 
Aguilar et al. 2014). The data suggest that models with high 
jjmax 10 TV or higher) are disfavored. We also note that 
the value j -py consistent with the naive estimate 

made from equating E^®® with 1 At this point it is clear 
that a pure OSNR scenario, which describes well the ^ GV 
- TV observations, cannot account for the CR hadronic flux 



Fig. 2.— Energy spectra of H (top) and He (bottom) multiplied by 
The solid lines indicate the model calculations. The contribution arising from 
OSNR (short-dashed lines) and from GSNR (long-dashed lines) are shown. 
The data are from PAMELA (Adriani et al. 2011), ATIC-2 (Panov et al. 
2009) and CREAM (Yoon et al. 2010). 


observed at ~ TV - PV rigidities. This is also the rigidity re¬ 
gion where the spectra are found to be harder. This considera¬ 
tion motivated us to introduce a second component for the CR 
hadron spectra at high energies, i.e., the GSNR component, 
representing the large-scale population of distant SNRs. Typ¬ 
ical parameters for GSNRs with strong shock and amplified 
magnetic fields are ui ^ 10 ®cms“^, B/kb 100p,G, and 
i?“®x ^ 5 PV. It is easy to see that, from these values, GSNRs 
are unable to produce and accelerate secondary or Li-Be- 
B. Eurthermore, the resulting CR spectra are totally insensi¬ 
tive to their exact values (and to the type of diffusion) so that 
the only relevant GSNR parameters are the source spectral in¬ 
dices. Eor both components, OSNR and GSNR, the slope a 
and theirnor normalization are chosen to match the data on 
primary spectra after propagation. The source parameter a 
is degenerated with the transport parameter 5, but the latter 
can been tested against the B/C ratio. As in Mertsch & Sarkar 
(2014) and related works, for Z = 1 we use a source spectral 
index steeper by 0.1 compared to that of heavier nuclei. This 
is a known issue, possibly ascribed to an A/Z-dependent in¬ 
jection efficiency in SNR shocks (Maikov et al. 2012). The 
relative source abundances are those adopted from previous 
studies (Tomassetti 2012; Tomassetti & Donato 2012) and we 
use the same values for the two SNR components. The con¬ 
tributions of the two components, determined from the data, 
are taken as 85 % for the OSNR and 15 % for the GSNR flux 
at 1 GeV/n, for all elements. Leptonic spectra from GSNR 
are expected not to contribute significantly to the high-energy 
flux, which is the case if these sources are placed at distances 
d >kpc (Delahaye et al. 2010). The data at >TeV energies 
require the GSNR spectra to be harder than those from the 
OSNR: we adopt an =4.1 and az>i =4.0. This is in fact 
encouraging, because the basic DSA predictions, supported 
by 7 -ray observations of young SNRs, favor a ^4.0 - 4.2 
(Blasi 2013). With this setup, in Eig. 2 we plot the model 
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Fig. 3.— Top: CR energy spectra of B and C multiplied by E^. Bottom: 
B/C ratio. The lines indicate the model calculations. In the top panel, the B 
components are for </)g + (short-dashed line), and (long- 

dashed line) and their sum (solid line). The C components are for (short- 
dashed line) and <j>^ (long-dashed line) and their sum (solid line). The data 
are from and PAMELA (Adriani et al. 2014), CREAM (Yoonetal. 2010), 
AMS-01 (Aguilar et al. 2010), and ATIC-2 (Panov et al. 2009). 

predictions for H and He fluxes as function of kinetic energy 
per nucleon. The two components of the flux are shown, i.e., 
split as (jiH = (t)H + iind cj^He = (I^He + ^He- The 
data are well described by the model, which interprets the 
TeV spectral hardening in terms of a smooth transition be¬ 
tween the OSNR and the GSNR component, but the fine struc¬ 
tures of the PAMELA data cannot be recovered. These sharp 
breaks seem also to be in contrast with the preliminary results 
of AMS. We eagerly await the AMS final results on H and 
He spectra at high energy, that will hopefully clarify how and 
where the spectral transition take places. 

We now come to the B/C ratio. Figure 3 shows the B and C 
fluxes (top) and the B/C ratio (bottom) compared with recent 
data, including the new data from PAMELA (Adriani et al. 
2014). The carbon flux, mostly primary, is also of the type 
</’C ~ </’c T ^C- experience a spectral hardening that is 

well reproduced by the model. The B spectrum, entirely sec¬ 
ondary, can be ideally split into cpB, = (j)% + </>]fwhere the 
first component is the one produced inside the OSNR, and the 
second arises in the ISM via collisions of heavier nuclei such 
as C, O, or Fe. Thus, the ISM component ^ can be split 

again into (produced by collisions of OSNR-emitted 

nuclei) and (by GSNR-emitted nuclei). In previ¬ 

ous OSNR-related works such as Cholis & Hooper (2014), 

the B/C ratio is always meant as (i/fg + (/'if 

markably, while the ratio (i/fg + starts rising 

at i? > lOOGeV per nucleon (as expected by a pure OSNR 
scenario), the ratio of the two-component model 4>b/4>C de¬ 
creases with energy in good agreement with the data. The 
trend of our B/C ratio is similar to the one from conventional 
propagation models (long-dashed line) where only GNRS 


component is considered. This effect can be understood from 
the top panel of the figure: at ^ lOOGeV/nucleon, when the 
ONSR component of B would become relevant enough to pro¬ 
vide a signature (i.e., 4% ^ </>b ^^°), the total fluxes of B and 
C are both dominated by the GSNR components, ' and 
4% respectively. Thus, in our two-component scenario, the 
B/C ratio does not constrain secondary production in SNR, as 
the presence of the GSRN component breaks the connection 
between the positron fraction and the secondary-to-primary 
nuclear ratios. We have tested the calculation using different 
values of But due to a degeneration with the OSNR 

parameter a, the parameter cannot be constrained much 
within the precision of the data. Our results for the B/C ra¬ 
tio seem to be quite robust: once accounting for the GSNR 
component, it decreases with energy as ~ E~^ (giving only 
little wiggles around the energy E k, see Fig. 3). 

For R^^^ ^ 1 TV or less, the spectra of H and He experience 
deviations from the power-law behavior below the TeV region 
that should be measurable by AMS. 

4. CONCLUSIONS 

Our work is motivated by two outstanding problems in 
CR physics: the origin of the rise in the positron fraction at 
E > lOGeV, and the nature of the ^TeV spectral hardening 
in CR hadrons. We revisited the ONSR scenario, proposed 
for the positron fraction anomaly, in order to account for 
the high-energy observations of CR hadron and nuclei. In 
OSNR models, secondary particles such as positrons and 
light-nuclei are produced and accelerated inside SNRs via 

hadronic interactions. In order to be able to match the 
data, these old-SNRs must have particular properties in 
terms of environmental parameters (such as strongly damped 
magnetization, or relatively high gas density) and predicts 
features in the B/C ratio that are not observed. In this Letter, 
we have argued that the OSNR scenario is incomplete for 
explaining the flux of CR hadrons at ^ TeV - PeV energies. 
The OSNR can account for the leptonic flux and for the 
GeV-TeV production of CR hadrons, but the flux at higher 
energies has to be provided by a population of distant and 
young sources that are able to efficiently accelerate CRs up 
to the knee. These sources are unable to produce secondary 
CRs (due to magnetic field amplification) and they do not 
contribute significantly to the high-energy leptonic flux 
(due to the large distance). Within this picture the spectral 
hardening of CR hadrons is interpreted as a signature of 
the transition between the local OSNR component and the 
Galactic ensemble. The spectra of all primary nuclei are 
predicted to harden. Taking into account the contribution 
of the two populations, we found that the predicted B/C 
ratio show no prominent signatures: it decreases with energy 
in accordance with the existing data. In conclusion, this 
generalized scenario may explains the absence of signatures 
in the B/C ratio while accounting for the observed signatures 
in primary CR hadrons. A quantitative inspection will 
be done with a proper modeling of leptonic and hadronic 
spectra arising from a realistic time-space distribution of their 
sources. It will be crucial, for achieving this goal, to have 
precision data on CR protons and nuclei in the energy region 
where the spectral transition takes place. 
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